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The ability of ionizing radiation to initiate genomic instability has been harnessed in the clinic where the localized 
delivery of controlled doses of radiation is used to induce cell death in tumor cells. Though very effective as a therapy, 
tumor relapse can occur in vivo and its appearance has been attributed to the radio-resistance of cells with stem cell- 
like features. The molecular mechanisms underlying these phenomena are unclear but there is evidence suggesting an 
inverse correlation between radiation-induced genomic instability and global hypomethylation. To further investigate 
the relationship between DNA hypomethylation, radiosensitivity and genomic stability in stem-like cells we have studied 
mouse embryonic stem cells containing differing levels of DNA methylation due to the presence or absence of DNA 
methyltransferases. Unexpectedly, we found that global levels of methylation do not determine radiosensitivity. In 
particular, radiation-induced delayed genomic instability was observed at the Hprt gene locus only in wild-type cells. 
Furthermore, absence of Dnmtl resulted in a 10-fold increase in de novo Hprt mutation rate, which was unaltered by 
radiation. Our data indicate that functional DNMTs are required for radiation-induced genomic instability, and that 
individual DNMTs play distinct roles in genome stability. We propose that DNMTS may contribute to the acquirement of 
radio-resistance in stem-like cells. 



Introduction 

Genome instability is a feature of cancer cells, manifesting in most 
cases as chromosomal instability (GIN), but in some endometrial 
and colorectal tumors as microsatellite or nucleotide-level insta- 
bility (NIN).^ This instability is an important part of oncogenesis 
and can be triggered by several environmental factors. Exposure 
to ionizing radiations is one such factor, which has been shown 
to promote malignancy in normal tissues.^ Radiation-induced 
genomic instability is characterized in vitro by several endpoints 
including elevated levels of non-clonal mutations, chromosomal 
aberrations and cell death occurring within a clonal cell popula- 
tion at a delayed period after irradiation.^ In vivo, the appearance 
of elevated levels of DNA damage, mutations and chromosomal 
instability can be detected in the un-irradiated progeny of irradi- 
ated animals. ^'^ 

Although the mechanism by which radiation-induced 
delayed genomic instability is propagated is poorly understood,^ 



several processes have been suggested as possible causes, includ- 
ing inflammatory type responses resulting in oxidative stress*^ 
and alterations in the balance of epigenetic factors such as DNA 
methylation.^'^ As a DNA-associated factor, DNA methylation 
is a promising candidate for the transmission and propagation 
of radiation-induced genomic instability.^ In vivo studies have 
demonstrated the occurrence of radiation-induced reductions 
in DNA methylation in the radiation target tissue, in a geno- 
type-, sex- and tissue-specific manner,^'^^ while also suggesting 
that the least methylated tissues tend to be more radiosensitive.^ 
In some instances, radiation-induced reductions in DNA meth- 
ylation have been coupled with reduced expression levels of the 
DNA methyltransferases,^^ and significant accumulation of DNA 
strand breaks. 

Aberrant hypomethylation, particularly when occurring at 
repetitive elements, including structural repeats and transpos- 
able sequences, is a common feature of cancer cells, and it is a 
major cause of genomic instability. Indeed, increased incidences 
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of recombination and rearrangements at sub-telomeric repeats 
and pericentromeric classical satellite sequences, ^^'^^ elevated lev- 
els of microsatellite instability,^^ reduced stability of triplet repeat 
sequences^^ and activation of endogenous retroviral elements,^^ 
have all been shown to be associated with hypomethylation. 
Furthermore, consistent with increased genomic instability in 
hypomethylated cells is the observation that mouse embryonic 
stem cells (mESCs) lacking the maintenance DNA methyltrans- 
ferase enzyme (DNMTl) display a 10-fold increase in mutation 
rate at Hprt and Tk gene loci.^^ Alongside global hypomethyl- 
ation, both tumor cells and DNMT3A knockout cells show local 
hypermethylation particularly at CpG islands. The mechanism 
for the hypermethylation is unknown, but apparently not caused 
by overexpression of DNMTs.^^ 

In the clinic, controlled delivery of ionizing radiation to 
tumor tissues is used to cause local DNA damage and subse- 
quent cell death making radiotherapy an extremely valuable 
and generally effective treatment for cancer.^^ However, sev- 
eral examples of radio-resistance and tumor relapse following 
radiotherapy have been reported, and it has been proposed that 
these may be due to the presence of a population of resistant 
cancer cells with stem cell features. It is unknown whether 
cancer stem cells can survive irradiation and/or if the radia- 
tion can induce stem cell characteristics in some of the cancer 
cells. Regardless of the mechanism, presence of radio-resistant 
cancer cells with stem cell characteristics, i.e., the ability to 
initiate a tumor, become evident at delayed times after the ini- 
tial irradiation. Cancer-initiating cells, or cancer stem cells, 
share several characteristics with embryonic stem cells, such as 
self-renewal capacities, the ability to form tumors, a common 
genetic program and similar DNA methylation profiles. ^^'^^ To 
further investigate the relationship between DNA hypometh- 
ylation, radiosensitivity and delayed genomic instability in 
stem cells, we have used an established panel of five mESCs 
with differing levels of DNA methylation due to the presence 
(wild-type line Jl) or absence of maintenance [Pnmtl'''Y^ 
or de novo (Dnrntda''' , Dnmtdh '' and Dnmt3a3b'^'Y^ DNA 
methyltransferases. For our studies we chose a radiation dose of 
3 Gy X-rays. This is a leukemogenic dose of ionizing radia- 
tion,^*^ sufficient to induce hypomethylation of murine bone 
marrow at a delayed time point post irradiation.^ Importantly 
for this study, 3 Gy is a biologically relevant dose, close to the 
fractions sizes routinely used in radiotherapy, and we reasoned 
that use of higher doses would have reduced the biological rel- 
evance of this study. 

Finally, analysis of cells immediately following radiation expo- 
sure was not included because the current study aimed to analyze 
instability at a delayed time point post radiation exposure rather 
than directly induced damage. The time points selected for the 
analyses were chosen to allow sufficient time for any changes 
in methylation or expression levels of genes to manifest. Time 
points were measured as population doublings (PDs) rather than 
days in culture to enable the comparison of the different mESC 
lines after the same number of cell divisions following irradiation 
with 0 or 3 Gy X-ray (see Fig. SI for the relationship between 
PDs and days in culture) . 



This study has allowed us to uncouple global DNA methyla- 
tion and radiation-sensitivity while revealing unique roles played 
by individual DNA methyltransferases in genomic stability. 

Results 

Global levels of DNA methylation in mESCs remain constant 
following irradiation. Global levels of DNA methylation in the 
panel of mESCs were assessed by HPLC-UV and found to be 
consistent with reported data.^^'^^ Relative to wild-type mESCs, 
functional absence of DNMTl substantially decreased genomic 
methylation, and while knockout of a single de novo methyl- 
transferase caused a very modest reduction, functional absence 
of both DNMT3A and DNMT3B led to a decrease in cyto- 
sine methylation that declined further with increasing passages 
(Fig. lA). Nevertheless, none of the mESCs displayed signifi- 
cant differences in cytosine methylation levels as measured at 
1 h, 24 h and 10 d post irradiation with 0 Gy (sham treatment) 
or 3 Gy X-ray (Fig. IB). There were no obvious abnormalities 
in the knockout mESCs with respect to growth rate, mor- 

phology or cell cycle distribution in comparison to the wild-type 
cell line, either before or after irradiation (Figs. SI and S2). The 
1 and 24 h time points were included in this analysis because a 
previous study of radiation-induced alterations in methylation 
levels reported that changes can occur as early as 6 h post irra- 
diation.^^ This time frame indicates that the causative mecha- 
nism is a driven response to radiation, rather than due to gradual 
loss of methylation through successive cell divisions. Thus, in 
analyzing whether radiation induced hypomethylation occurred 
in the 5 mESC lines, we included the early time points to assess 
whether any alterations observed were due to passive loss during 
cell division or by driven demethylation. The results, however, 
indicated that no alterations in methylation level occurred in 
the mESCs post irradiation. The inclusion of early time points 
in other experiments was thus deemed unnecessary as the main 
focus of this study were the delayed effects of radiation exposure. 

Radiosensitivity in mESCs is independent of global levels 
of DNA methylation. Long-term clonogenic survival assays on 
mESCs after irradiation with X-ray doses increasing from 0-7 
Gy revealed a dose-dependent decrease of the surviving fraction 
(SF) in all five lines (Fig. IC). Similarly to the data reported 
by other groups working with mESCs, none of the cell lines 
showed any evidence of a shoulder, but instead displayed a near 
exponential survival curve (Fig. IC), as confirmed by the fact 
that the (3 values obtained by fitting a linear quadratic model to 
the data are all very close to zero (Fig. ID). Analysis of the sur- 
vival curves showed no significant difference between the wild- 
type cell line and any of the hypomethylated D;^;/?^ knockout cell 
lines. However, the Dnmt3a3h''' cell line was significantly more 
radiosensitive than Dnmt3a''' and Dnmt3h''' (SF3Gy, p = 0.005 
and p = 0.015, respectively). It must be noted that Dnmt3a3h''' 
cells were analyzed at passage 20-25, when these cells have simi- 
lar levels of DNA methylation to Dnmtl '' cells (Fig. lA). As the 
two most hypomethylated cell lines (Pnmt3a3h''' and Dnmtl ' ) 
display similar clonogenic survival responses to radiation as 
the wild-type cell line, it seems that hypomethylation alone is 
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Figure 1. Relationship between DNA methylation and irradiation. Genomic DNA cytosine methylation measured by HPLC-UV of un-irradiated mESC 
samples and mouse control tissues before (A) and after irradiation (B). Columns represent the average percentage of cytosine bases that were methyl- 
ated. Error bars represent the standard error of the mean (SEM). Results were combined from multiple independent experiments. The number of 
replicate measurements for each sample is displayed above the error bar. (C) Dose-response curve of the ESC lines constructed using the long-term 
clonogenic assay after irradiation with 0-7 Gy X-rays. Each data point represents the mean, and error bars represent the SEM, of at least three indepen- 
dent experiments. (D) Linear quadratic parameters a and p were calculated by fitting the data displayed in (C) with a linear quadratic equation using 
SPSS and are shown in columns 3 and 5. Column 2 indicates the plating efficiency for each clone and column 4 and 6 the standard errors (SE) for a and 
P respectively. The survival fraction at 2Gy (SF2) inferred from the data are shown in the last column. 



not sufficient to increase radiosensitivity. Thus, the DNMT3 
family of enzymes may contribute to the determination of 
radiosensitivity in mESCs, and the absence of a single de novo 
methyltransferase, particularly DNMT3A, may result in slight 
radioresistance. 

Delayed effects of radiation as measured by DNA damage 
or oxidative stress are comparable in v^^ild-type and hypometh- 
ylated mESCs. Radiation induced delayed genomic instability 
is characterized by the persistent expression of elevated levels 
of non-clonal mutations or aberrations at a delayed period after 
irradiation, which are detected in the un-irradiated progeny of 
the irradiated cells^^ or animals^ several generations after irradia- 
tion.^ The persistent presence of elevated levels of DNA damage, 
manifesting as strand breaks and oxidative damage in the un- 
irradiated progeny of irradiated cells, is an accepted indicator of 
ongoing genomic instability, and the Comet assay has success- 
fully been used to detect such ongoing instability in multiple 



previous studies. ^'^^ We have used the Comet assay as another 
measure of delayed instability and we thus looked at late time 
points to analyze those cells surviving, and "recovered from," 
the effect of irradiation, i.e., the un-irradiated progeny of the 
irradiated cells. 

Average levels of single-strand breaks (SSBs), double-strand 
breaks (DSBs) and alkali-labile sites (ALS) were determined in 
wild-type, Dnrntl'^' and Dnmt3a3b'^' mESCs at 23-25 popula- 
tion doublings (PDs) after treatment with 0 or 3 Gy X-irradiation 
using the alkaline comet assay. The average level of endoge- 
nous and/or induced DNA damage was found to be very simi- 
lar between all of the cell lines and treatment conditions (p > 
0.07) (Fig. 2A). The Fpg-modified comet assay was then used to 
assess the level of oxidative damage in wild-type, Dnmtl '' and 
Dnmt3cL3h''' mESCs. After correction for damage levels observed 
in the non-enzyme treatment controls, there was no significant 
difference in the average level of FPG-sensitive lesions between 
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Figure 2. Analysis of single strand, double strand and alkali labile site 
DNA damage (A) or damage resulting from oxidative stress (B) in wild- 
type, Dnmtr^-. DnmtSaSb-^- mESC lines, 23-25PDs post 3 Gy X-irradiation 
or sham treatment as measured by percentage of tail DNA. Results are 
expressed as an average of three independent experiments, each com- 
prising 200 separate cell measurements. Error bars represent the SEM of 
the averages of the three independent experiments. 



any of the cell lines (p > 0.07), or between the 3 Gy and 0 Gy 
(sham) treatment conditions (p > 0.06), indicating no differences 
in the level of damage induced by long-term oxidative stress (Fig. 
2B). 

Radiation-induced delayed chromosomal instability is inde- 
pendent of DNA methylation levels. Next, the frequencies of 
cytogenetic aberrations were assessed in wild-type, Dnrntl'^' and 
Dnrntdadh '' mESC lines and the structural cytogenetic aberra- 
tions observed were grouped according to defined classifications^^ 
at 10-14 (Fig. 3A) and 23-25 (Fig. 3B) population doublings 
post irradiation with 0 Gy or 3 Gy X-rays (see Table 1 for 
details). A mid-time point of 10-14 PDs was used in the cytoge- 
netic analyses, in order to see whether levels of aberrations varied 
over time post irradiation and the results of cells examined post 
exposure were always compared with the results for un-exposed 
control cells for the same time point. There were no significant 



differences in the frequency of total structural aberrations or of 
any of the individual aberration types between any of the cell 
lines, treatment groups or time points, as calculated using the 
chi-square test (p > 0.14 and p > 0.11, respectively). The average 
number of aberrations per metaphase was also not significantly 
different between any of the cell lines, treatment conditions or 
time points (p > 0.12, data not shown). 

A high ratio of chromatid: chromosome type aberrations 
is accepted as indicative of ongoing radiation-induced chro- 
mosomal instability in cultures of cells. ^^'^^ In this study, the 
chromatid: chromosome type aberrations ratio observed in the 
wild-type cell line was comparable at both time points (Table 1). 
In addition, this ratio did not differ significantly between the 
3 Gy and sham treatment groups for any of the mESC lines at 
both time points. 

Possible relationship between methylation levels and aber- 
rations involving the satellite arms. An aberration involving 
the appearance of abnormally long satellite arms was observed 
almost exclusively in the Dnmtl '' cell line at the 23-25 PD time 
point (Fig. 3C). It occurred once in the 3 Gy-treated wild-type 
cells, eight times in the 3 Gy-treated Dnmtl '' cells, and five times 
in the sham-treated Dnmtl '' cells (see Table 1 for frequency of 
abnormalities). As it could not be determined whether these aber- 
rations were clonal, a conservative approach was taken and mul- 
tiple appearances within the same cell population were counted 
as a single aberration. In addition, several complex exchanges 
involving the satellite arms of three or more chromosomes were 
observed solely in the hypomethylated Dnmtl '' and Dnmt3a3b''' 
mESC lines, at both time points and irrespective of the treat- 
ment condition (Fig. 3D). Assuming that the long satellite arm 
aberration was clonal, the difference in total frequency of aber- 
rations involving associations or exchanges of the satellite arms 
in comparison to wild-type cells was not statistically significant 
for Dnmtl '' mESCs (p = 0.17), while it displayed a significant 
increase in Dnmt3a3b''' cells (p = 0.006). 

DNMTs are required for radiation-induced delayed genomic 
instability. As we did not detect any highly significant differ- 
ences in delayed instability cytogenetically between the mESC 
lines, we performed the analysis at a higher resolution level by 
comparing mutation rates at a single gene locus, Hprt. Locus spe- 
cific delayed genomic instability was assessed in mESCs cultured 
in the presence of 6TG and the mutation rate at the Hprt locus 
was measured 23-25 population doublings after irradiation with 
0 Gy and 3 Gy X-rays (sham and 3 Gy in Fig. 4A, respectively) . 
This permitted analysis of the de novo mutation rate of each cell 
line in addition to the radiation-induced delayed mutation rate 
(see Table 2 for details of selection). The appearance of individ- 
ual colonies in the 6 TG selection dishes revealed that no founder 
mutations were present in any of the clonal populations. 

Mutation rates were calculated using Luria-Delbriick fluc- 
tuation analysis. The de novo mutation rates of the Dnmt3a''', 
Dnmt3b''' and Dnmt3a3b''' mESCs were not significantly differ- 
ent from the wild-type mESC line (p > 0.2) (Fig. 4A). Therefore, 
global hypomethylation does not appear to influence mutation 
rate at the Hprt locus in mESCs. However, in agreement with pre- 
vious reports,^^ the Dnmtl '' mESCs displayed a 10-fold elevated 
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Figure 3. Classifications of structural cytogenetic aberrations as defined by Savage. The results are displayed as the number of each class of aberra- 
tion observed per 100 metaphases. Graph (A) displays the results 10-14 PDs post treatment; graph (B) displays the results 23-25 PDs post treatment. 
(C) Mitotic spreads of Giemsa-stained chromosomes from the DnmtV^' ESC line 23-25 PDs after 0 Gy treatment, showing a chromosome with very 
long satellite arms. Immediately after exposure to 3 Gy X-rays or sham treatment, cells were seeded into the clonogenic assay. Surviving colonies 
were counted 12-14 d after seeding. (D) Dicentric, tri-centric and tetra-centric chromosome-type aberrations observed in Giemsa-stained metaphase 
spreads of Dnmt7"/"and Dnmt3a3b'^' mESCs 23-25PDs post OGy (sham) or 3 Gy treatment. 



de novo mutation rate in comparison to the wild-type mESCs (p 
= 2.7 X lO'*^), indicating the importance of DNMTl for genomic 
stabihty in mESCs (Fig. 4A). Unexpectedly, the de novo muta- 
tion rate of Dnrntl'^' mESCs was significantly higher than that 
of Dnmt3a3h-'- mESCs (p = 3.2 x 10'^). Dnmtdadh '- mESCs 
were at passage 20-26 at the time of seeding into selective media 
and therefore had a very similar low global methylation level 
to Dnmtl '' mESCs (Fig. 1). Thus, it appears that the elevated 
mutation rate detected in the Dnmtl '' mESC line is not a result 
of global hypomethylation but a result of the absence of a prop- 
erty or function specific to DNMTl or of loss of methylation at 
genomic regions specifically targeted by DNMTl. 

Wild- type mESCs displayed a significantly elevated muta- 
tion rate in the 3 Gy X-irradiated clones in comparison to the 
sham treated control clones at 23-25 population doublings 
post treatment (p = 1.6 x 10'^). This indicates the occurrence 
of radiation-induced delayed genomic instability. However, none 



of the Dnmt knockout mESCs showed a significantly elevated 
mutation rate in the 3 Gy X-irradiated clones in comparison to 
their corresponding sham treated controls (p > 0.4) (Fig. 4A). 
Thus, Dnmt knockout mESCs showed no evidence of radiation- 
induced delayed genomic instability. These results indicate that 
absence of the DNA methyltransferase enzymes, or disruption 
of the normal methylation pattern, may potentially inhibit the 
mechanism behind radiation-induced delayed genomic instabil- 
ity at the nucleotide level. 

Increased retrotransposition in Dnmtl ' mES cells. To inves- 
tigate the mechanisms that might be responsible for the differing 
mutation rates between the ESC lines, the types of mutations 
present in the Hprt gene were characterized using exonic PCRs 
on DNA from 6-TG resistant clones from each cell line and treat- 
ment condition. The spectra of mutations detected in each cell 
line after sham or 3 Gy treatment are displayed in Figure 4B; 
Figure 4C displays the numbers of each mutation type identified 



896 



Epigenetics 



Volume 7 Issue 8 



Table 1. Absolute numbers of each cytogenetic aberration type observed at 10-14 PD and 23-25 PD post treatment witli OGy or 3Gy X-rays 













Number of Chromosome-type Aberrations (Frequency) 


Number of Chromatid-type Aberrations (Frequency) 


■ Other 












Simple Interchange 


Simple Intrachange 


Simple Interchange 


Simple Intrachange 


Population 
doubling 


Cell line/condition 


No. of 

Metaphases 

Scored 


No. of 
aberrant 
metaphases 
(frequency) 


No. of 

aberrations 

(frequency) 






CR/ ins del 


Discontinuity 






CR/ ins del 


Discontinuity 


No. of 

Compl6x 

Satellite 

/ Exchanges 
(Frequency) 


Ratio of 
Chromatid: 

Type 

Aberrations 




Wild Type Sham 


119 


15(0.126) 
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0 


5 (0.042) 


0 


0 


0 
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Wild Type 3Gy 
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21 (0.088) 


23 (0.097) 


0 


3 (0.013) 


0 


8 (0.034) 


0 


0 


0 


12 (0.050) 


0 


1.1 


10-14PDS 


Dnmt1-/- Sham 


118 


16(0.136) 


22 (0.186) 


0 


2 (0.017) 


0 


13 (0.110) 


0 


0 


0 


6 (0.051) 


1 (0.008) 


0.4* 


Dnmtl-/- 3Gy 


248 


30 (0.121) 


36 (0.145) 
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0 
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0 


0 
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Dnmtl-/- 3Gy 
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1 (0.008) 


0 
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1.1 
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119 
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0 
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1 (0.008) 


0.8 
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One liundred to one liundred and twenty metaphases were scored for each cell line in each treatment condition. *The ratio of chromatid:chromosome 
type aberrations is significantly different between the sham and 3 Gy treatment groups at 10-14PDs for Dnmtl-/- cells (X^ p=0.024 after Bonferroni cor- 
rection). RT, reciprocal translocations; Die, dicentric; CR, centric ring; insdel, insertion or deletion. 



in sham treated and 3 Gy irradiated cells. The details of mutants 
obtained and analyzed by PGR are summarized in Table SI. 
Multiple deletions of non-contiguous exons were observed in 
sham treated clones from the two most hypomethylated mESG 
lines {Dnrntl'^' and Dnmt3a3h''') and in clones derived from all 
of the mESG lines after irradiation with 3 Gy X-ray (Fig. S3). 
Hypomethylated Dnmtl '' mESGs, which had the highest de 
novo mutation rate, also displayed the widest variety of muta- 
tion types, indicating that several different mutational processes 
may occur in this cell line. Dnmtl '' cells were the only line to 
demonstrate the occurrence of active retrotransposition, indi- 
cated by insertions involving two different subfamilies of non- 
autonomous SINE retrotransposons (B2 and Bl/Alu) 5' of exon 
3 and exon 6 respectively (Fig. S4). 

It is possible that active retrotransposition was not observed 
in the other lines due to their lower spontaneous mutation rates. 
However, the high de novo mutation rate of Dnmtl '' mESGs 
was further demonstrated by the occurrence of second mutation 
events during expansion of the clones. For example, one of the 
many colonies carrying the Bl SINE insertion was also found 
to have a deletion of exon 4, which must have arisen at a delayed 
period after the insertion. As SINEs require the LINE-1 encoded 
ORF2 protein for their retrotransposition,^*^ LINE-1 methylation 
status was tested in these cells. Although all the mESGs are hypo- 
methylated relative to control tissues, Dnmtl '' and Dnmt3a3b''' 
LINE-ls were slightly less methylated than wild-type mESGs 
(Fig. S5), suggesting that LINE-1 hypomethylation per se is not 
sufficient to promote retrotransposition as Dnmtd^il 3h''' mESGs 
do not show such an effect. 

Discussion 

The work described in this study was aimed to examine the 
relationship between DNA methylation, radiosensitivity and 
delayed genomic instability in mouse embryonic stem cells. 
Rather than using inhibitors of DNA methyltransferases such as 



5-azacytidine, which have a broad range of effects, we addressed 
the roles of maintenance (DNMTl) and de novo (DNMT3A 
and DNMT3B) DNA methyltransferases using a panel of mEGS 
that is well established for the studies of activities of individual 
DNMTs. This panel consists of five mESG lines including the 
wild-type line JP^ and four of its derivatives (Dnmtl''', Dnmt3a''', 
Dnmt3b''' and Dnmt3a3b''') in which the DNA methyltransfer- 
ases were inactivated by targeted deletions within their catalytic 
domains. ^^'^^ 

There are reports describing alterations in methylation levels 
post 3Gy X-irradiation in somatic tissues in vivo,^'^^ and a cor- 
relation between global DNA hypomethylation and increased 
radiosensitivity typical of somatic cell lines. ^^'^^ However, radi- 
ation-induced methylation was observed to vary in a strain-, 
tissue- and sex-specific manner. We have now shown that in 
mESGs, exposure to 3 Gy X-irradiation does not lead to changes 
in DNA methylation, and that radiosensitivity is independent 
of DNA methylation levels. These findings may reflect unique 
properties of mESGs in response to radiation that distinguish 
them from somatic cells, though further studies will be required 
to confirm this. Given the cell- and tissue-type specificity of these 
phenomena,^^ our finding may represent a property of mESGs 
reflecting their hypomethylation in comparison to most somatic 
cells. Gonsistent with previous observations, all five mESG lines 
displayed nearly exponential survival curves,^^ which may reflect 
their increased propensity to undergo apoptosis in response to 
DNA damage.^^ Yet, the de novo methyltransferases, DNMT3A 
and DNMT3B, may play a role in modulating sensitivity to 
X-rays in mESGs, as their absence seems to have a modest radio- 
protective effect. 

Although levels of radiation-induced delayed chromosomal 
instability were comparable in all the mESGs lines before and 
after irradiation, aberrations and multicentric exchanges /associa- 
tions involving the satellite arms of chromosomes were present 
at an increased frequency in hypomethylated mESGs {Dnmtl''' 
and Dnmt3a3b''') compared with wild-type cells (Table 1). It is 
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Table 2. Details of the selections carried out to determine the mutation rate at the Hprt gene for each cell line and treatment condition 



Cell line and 
condition 


Number of independent 
selections 


Number of cells screened 
per selection 


Average plating 
efficiency 


Number of mutations per 
selection^ 


Mutation rate"^ 


Wild type sham 


25 


4 X 10^ 


55% 


0(22), 1(2), 3(1) 


5.7 X 10"^ 


Dnmt17-sham 


19 


4x 10^ 


52% 


0(6), 1(7), 2(2), 3(1), 4(1), 
12(1), 97(1) 


4.4 X 10-7 


Dnmt3a3b 
--/-- sham 


25 


4x 10^ 


51% 


0(22), 1(2), 2(1) 


6.0 X 10-^ 


1 Dnmt3a7-sham 


19 


4x 10^ 


38% 


0(16), 1(2), 2(1) 


1.1 X 10-7 


DnmtSbV- sham 


22 


4x 10^ 


51% 


0(16), 1(1), 1(4), 16(1) 


8.7 X 10-« 


Wild Type 3 Gy 


20 


4x 10^ 


49% 


0(11), 1(5), 2(1), 3(1), 9(1), 
13(1) 


2.8 X 10-7 


Dnmt1-/-3 Gy 


19 


4x 10^ 


49% 


0(7), 1(5), 2(2), 3(2), 6(2), 
29(1) 


4.5 X 10-7 


Dnmt3a3b 
--/-- 3 Gy 


20 


4x 10^ 


45% 


0(16), 1(3), 6(1) 


1.2 X 10-7 


Dnmt3a-/-3Gy 


20 


4x 10^ 


41% 


0(17), 1(2), 2(1) 


9.5 X 10-^ 


Dnmt3b-/- 3 Gy 


21 


4x 10^ 


48% 


0(18), 2(1), 7(1), 8(1) 


8.5 X 10-« 



+The number in brackets indicates the number of selections in which the adjacent number of drug-resistant colonies was observed. ^Calculated using 
Luria Delbruck fluctuation analysis and corrected for plating efficiency. 



possible that increased frequencies of rearrangements or homolo- 
gous recombination at the major and minor satellite repeats occur 
due to hypomethylation. Repetitive sequences, such as retroviral 
elements and satellite repeats, contain the majority of the methyl- 
cytosine in mammalian cells,^^ can be hotspots for recombina- 
tion^^ and are hypomethylated in patients carrying DNMTl 
mutations .^^ Gonzalo et al. demonstrated that Dnmtl '' and 
Dnmt3cL3h''' mESCs display greater heterogeneity of telomere 
length and significantly increased frequencies of telomere sister 
chromatid exchange (T-SCE). Interestingly, exchanges between 
hypomethylated satellite sequences, resulting in the formation of 
multi-radial chromosomes, are frequently observed in lympho- 
cytes from patients with ICF syndrome, which carry mutations 
in DNMT3B}^ Further work is required to determine whether 
the frequency of recombination throughout the genome is ele- 
vated in the hypomethylated Dnmtl '' and Dnmt3a3h''' mESCs. 

Sham treated Dnmtl '' mESCs displayed an approximately 
10 -fold elevation of de novo mutation rate at the X-linked Hprt 
gene locus in comparison to that of wild-type mESCs, while all 
the other cell lines did not differ from wild- type cells. A previous 
report has described an increase in mutation rates in Dnmtl '' 
cells; however, only Dnmtl '' cells were investigated, without any 
irradiation treatment, and it was concluded that the increased 
mutation rate was a result of DNA hypomethylation.^^ We have 
now demonstrated that this elevated mutation rate is not simply a 
result of global hypomethylation. If that were the case, a compa- 
rable increase in de novo mutation rate would be expected for the 
similarly hypomethylated Dnmt3cL3h''' mESC line. Instead, the 
mutation rate in this cell line does not differ from that of wild- 
type cells. Rather, the fact that the increase in mutation rate is 
specific to Dnmtl '' cells suggests that it is caused by the absence 
of methylation in regions specifically targeted by DNMTl or of 
other properties or functions specific to DNMTl. Such functions 



may include the various roles played by DNMTl in DNA rep- 
lication and repair. DNMTl, but not DNMT3A or DNMT3B, 
is recruited by PCNA to the replication fork during S phase of 
the cell cycle, and also to sites of DNA repair.^^ Furthermore, 
DNMTl is thought to play a role in strand discrimination dur- 
ing mismatch repair (MMR).^^ Thus, it is possible that absence 
of functional DNMTl or the lack of methylation of genes and 
regions targeted by DNMTl may hinder the efficient function- 
ing of DNA repair, causing an increase in mutation rate. In 
addition, mobilization of retrotransposons in Dnmtl '' cells may 
contribute to the overall increase in mutation rate, an observation 
consistent with recent evidence that degradation of DNMTl fol- 
lowing exposure to benza(a)pyrene leads to LI reactivation.^*^ A 
hypothesis that remains to be explored is whether there is also an 
elevated mutation rate in patients carrying DNMTl mutations 
that lead to late onset, neurological conditions.^^'^^ 

Non-contiguous deletions at the Hprt locus were seen in the 
most hypomethylated lines (DnmtT'' and Dnmt3cL3h''') before 
irradiation, and in all lines post irradiation. We propose that 
their appearance may reflect increased homologous recombina- 
tion (HR). Ionizing radiation is known to induce a persistent 
increase in homologous recombination,^ which can be detected 
as sister chromatid exchanges (SCEs).^^ Furthermore, Rugo et al. 
(2011) found that naive mESCs exhibited a significant increase 
in the frequency of SCEs when they shared media with descen- 
dants of irradiated mESCs, and that SCEs were increased in cells 
that shared media with un-irradiated triple DNMT knockout 
mESCs, compared with cells that shared media with un-irradi- 
ated wild-type mESCs. Further studies will be required for the 
detection of SCEs in our system. 

In contrast to what we observed in wild-type mESCs, none 
of the Dnmt knockout cells displayed radiation-induced delayed 
genomic instability, a direct demonstration that specific DNA 
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0 


0 
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1 


2 
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0 


2 


6 
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0 
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0 


0 


0 


0 


2 
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0 


0 
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0 


1 
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2 
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Figure 4. (A) Hprt gene mutation rates 23-25 population doublings after 3Gy X-irradiation or 0 Gy (sham) treatment. Columns represent the muta- 
tion rate, calculated from 19-25 individual clonal expansions using Luria-Delbruck fluctuation analysis and corrected for plating efficiency. Error bars 
represent the 95% confidence intervals. DnmtSaSb'^' mESCs were p#20-26 at the time of selection for mutants. (B) Spectrum of functional Hprt gene 
mutations observed in mESC lines 23-25 population doublings after 3 Gy X-irradiation or 0 Gy (sham) treatment. (C) Tables detailing the mutations 
identified by exonic PGR that were used to construct the pie charts in (B). PCRs were performed using primer pairs designed to amplify each exon (1-9) 
of the Hprt gene and exon 2 of the unlinked K-ras gene as a control. When a mutation was observed in multiple colonies arising in the same clonal 
population it was counted only once. "Other mutations" include mutations that could not be detected using the PCR-based method, such as point 
mutations, frame shift mutations and epigenetic alterations. 



methylation patterns and/or the DNA methyltransferases are 
required for this phenomenon to occur. These findings are largely 
congruent with recently published data obtained studying the 
bystander effect,^^ and indicate that active DNMTs are required 
for persistent radiation-induced genomic instability. It has been 
proposed that alterations in DNMT activity or the DNA meth- 
ylation profile could potentially alter the secretion profiles of cells 
for factors such as TNFa, NO and TGF|3, which impact genomic 
stability. The inheritance of such altered methylation patterns 
could potentially result in propagation of the instability phe- 
notype through successive generations (as seen in un-irradiated 
Dnrntl'^' mESCs), or lack of appearance of the instability phe- 
notype after radiation exposure (as seen in all Dnmt knockout 
mESCs exposed to 3 Gy X-rays). Further investigation of the 
expression levels of such factors in the mESCs would be neces- 
sary to determine which, if any, of these possible scenarios is true. 
The data point to an apparent paradoxical nature of DNMTl 
loss, which can both promote (un-irradiated Dnmtl '' mESCs) 



and prevent (irradiated Dnmtl '' mESCs) delayed genomic insta- 
bility. However, the lack of increase in genomic instability levels 
in Dnmtl '' cells following irradiation may be a reflection of the 
already high mutation rates in un-irradiated cells. Future experi- 
ments, including the rescue of knockout cells with the respec- 
tive Dnmts and treatments with demethylating agents and/ 
or specific DNMT inhibitors, will be required to investigate 
whether absence of functional DNMTs or disruption of the nor- 
mal methylation pattern inhibit the causative mechanism behind 
radiation-induced delayed genomic instability or the appearance 
of the effects. 

In summary, we have provided experimental evidence sug- 
gesting that radiosensitivity or delayed genomic instability in 
mESCs is independent of global levels of genomic DNA hypo- 
methylation, but that individual DNMTs play distinct roles in 
these responses, whereby the de novo DNA methyltransferases 
may influence radiosensitivity and are required for radiation- 
induced genomic instability, and functional DNMTl is required 
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to maintain the stability of genes such as Hprt. As mouse embry- 
onic stem cells share several characteristics with cancer stem 
cells, our data support the hypothesis that radioresistance and 
genomic instability in cancer may in some instances be attributed 
to changes in activity of DNMTs (e.g., ref 49). Analysis of these 
factors in normal and tumor tissues will further our understand- 
ing of the conditions /processes required for the manifestation 
of radiation induced genomic instability, which may eventually 
lead to oncogenesis and/or radioresistance. At the same time, this 
experimental system will allow a further characterization of the 
mechanisms by which lack of DNMT activities contributes to 
these phenomena and the identification of methods to modulate 
radiosensitivity in tumors and normal tissues, thus contributing 
to the design and optimization of future therapies. 

Materials and Methods 

Sequences of oligonucleotides can be found in Tables S2 and S3. 

mESCs culture and irradiation. mESCs cell lines were main- 
tained on gelatinized plates in DMEM high glucose (Invitrogen) 
with ESGRO (Millipore) as previously described.^^ mESCs were 
irradiated with 1-, 3-, 5-, 7- or 10-Gy of X-rays delivered at 0.5 
Gy min'^ (250 kV constant potential, HLV 1.5 mm Cu), using a 
Pantak industrial X-ray machine, and media changed immedi- 
ately after irradiation. 

DNA methylation analysis by HPLC-UV. Genomic DNA 
was extracted using QIAamp DNA Mini Kit (Qiagen) follow- 
ing treatment of cell pellets with RNase A (Roche Diagnostics, 
lOmg/ml) and RNase Tl (Roche Diagnostics, 100,000 U/ml). 
Global levels of DNA cytosine methylation were characterized 
using high performance liquid chromatography with UV detec- 
tion (HPLC-UV) as previously described.^^ 

Clonogenic assays. Clonogenic assays were performed as pre- 
viously described^^ by seeding 10^, 2 x 10^, 4 x 10^ and 5 x 10^ 
cells for 0, 1, 3, 5 and 7 Gy treated mESCs respectively. Plating 
efficiency (PE) was assessed by the proportion of cells seeded into 
the OGy well that survived to produce colonies. Statistical analy- 
sis was performed using a two-tailed T-test. 

Comet assays. Alkaline Comet Assays were based on previ- 
ously described protocols^^ and performed 23-25 population 
doublings post 3 Gy X-irradiation or sham treatment. Images 
were analyzed with the Komet Analysis software version 5.5 
(Andor Bioimaging). Two hundred cells were analyzed per sam- 
ple per experiment. Percentages of the cellular DNA within each 
comet tail were calculated (% tail DNA), and results expressed as 
the mean and standard error of the mean (SEM) of three inde- 
pendent experiments. Control slides prepared using wild-type 
mESCs irradiated with 0, 5 and 10 Gy X-Rays immediately prior 
to the Comet assay were included to ensure consistency between 
experiments. 

The relative levels of oxidative purine base damage were moni- 
tored using the human formamidopyrimidine DNA glycosylase 
(fpg) comet assay (fpg comet) as described previously for the 



human 8-oxoguanine DNA glycosylase 1 (hOGGl) comet assay 
(hOGGl comet) with the following modifications: after lysis 
the slides were washed once with distilled water and immersed 
in three changes of enzyme digestion buffer [40 mM HEPES, 
0.1 M KCl, 0.5 mM EDTA and 0.2 mg/ml bovine serum albu- 
min (pH 8.0)], for 5 min each time, on ice. Fpg (Sigma-Aldrich) 
was added to the gel (50 |xl/gel) at 1/500 dilution; gels were 
covered with a coverslip and incubated in a humidified chamber 
at 37°C for 30 min; the coverslips were removed and the slides 
were placed in a horizontal electrophoresis tank. From this step 
onwards, the assay was performed as described.^^ DNA damage 
was expressed as the percentage of DNA in the comet tails. 

Cytogenetic aberrations. Cells were treated with 0.03 |JLg/ml 
KaryoMAX colcemid solution (GIBCO) eight hours after irra- 
diation with 0 or 3 Gy X-ray and slides prepared for cytogenetic 
analysis as previously described.^^ Scoring was performed with 
bright field function and a green filter on a Carl Zeiss Microscope 
(Axioskop 2). Pictures were taken using Ikaros (MetaSystems) 
or Smart Capture X (Digital Scientific UK) imaging software. 
P values were calculated using the chi-square test. 

Hprt mutation analysis. Analysis of mutation rate at Hprt on 
0 Gy or 3 Gy X-rays treated cells was performed according to the 
6-thioguanine (6TG) selection protocol as described.^^ Mutation 
rates were calculated using the maximum likelihood method of 
Luria-Delbriick fluctuation analysis implemented in Salvador 
v2.3.^^ Exonic mutations were studied as previously described.^*^ 
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